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Abstract 

The supersymmetric extension of the standard model suffers from a problem of 
baryon-number violation. Discrete (and global) symmetries introduced to pro- 
tect the proton are unstable under gravitational effects. We add a gauged U(l)x 
to the standard model gauge group Gsm and require it to be anomaly-free. As 
new (chiral) superfields we only allow GsM-singlets in order to maintain the 
good unification predictions. We find the most general set of solutions for the ra- 
tional singlet charges. We embed our models in local supersymmetry and study 
the breaking of supersymmetry and U(l)x to determine Mx- We determine 
the full non-renormalizable and gauge invariant Lagrangian for the different so- 
lutions. We expect any effective theory to contain baryon- and lepton-number 
violating terms of dimension four suppressed by powers of Mx /Mpi. The power 
is predicted by the U(l)x charges. We find consistency with the experimental 
bounds on the proton lifetime and on the neutrino masses. We also expect all 
supersymmetric models to have an unstable but longlived lightest supersym- 
metric particle. Consistency with underground experiments on upward going 
muons leads to stricter constraints than the proton decay experiments. These 
are barely satisfied. 
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1 Introduction 

When incorporating supersymmetry into the Standard Model one immediately runs 
into a problem. The Standard Model conserves baryon- (B) and lepton-number (L) 
automatically and higher dimensional AB, AL ^ operators are suppressed by at 
least 4 powers of the scale of baryon- or lepton-number violation. However, in su- 
persymmetry the most general interactions involving the Standard Model (super-) 
fields and invariant under the Standard Model gauge group 

Gsm = SU(3) C <8> SU(2) L ® U(1) Y , (1) 
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include AB, AL ^ terms of dimension 4 J2], |, J| 



(LLE) F , (LQD) F , (UDD) F , {LH) f . (2) 

Here L, Q are the left-handed lepton and quark superfields respectively.^] E, U, D are 
the corresponding right-handed superfields. Their Gsm quantum numbers are given 
below in Eq.(|6]). There are also dim 5 terms which can lead to dangerous levels of 
proton decay [I], |2j 

(LLHH) F , (QQQL) F , (UUDE) F , (LHHH) F , 

(QQQH) F , (HQUE) F (QUL*) D , (UD*E) D , (3) 

(HHE*) D , (H*HE) D . 

With the gauge structure restricted to that of the minimal supersymmetric stan- 
dard model (MSSM) the LH term in Eq.(0) can be rotated away through a field 
redefinition. It is then absorbed in the LLE and LQD terms as well as those of the 
MSSM. However, if extra gauge symmetries are present which distinguish between 
L and H this is no longer true. This gauge symmetry could prohibit the first two 
terms of Eq.(|2]) as well as the fiHH term of the MSSM while allowing a kLH term. 
This can then no longer be rotated away. In the low-energy effective theory this extra 
symmetry is broken. If it is broken at sufficiently high energy no effects are observable 
and the LH term can again be rotated away. We thus retain the LH term in our 
discussion and decide if there is a remnant effective LH term. 

The dim 4-terms together lead to an unacceptable level of proton decay. Thus the 
symmetry of the SM must be extended to 

G = G S m®G, (4) 

such that G guarantees the longevity of the proton. Discrete, global, and gauge 
symmetries have been considered. The first solution to the problem of proton decay 
was to introduce the discrete symmetry R-parity [||, [3], Q which prohibits all the terms 
in (Q). Later other discrete symmetries were considered, in particular baryon-parity 
|| which only prohibits the term UDD. This is sufficient to protect the proton. 
However, Krauss and Wilczek later showed that discrete symmetries are violated 
by gravitational effects unless they are remnants of a broken gauge symmetry. So only 
very special discrete symmetries, discrete gauge symmteries, are viable candidates for 
G. Ibanez and Ross || |6|] translated the requirement of an anomaly-free gauge theory 
into a set of conditions on the remnant discrete symmetry. They then systematically 
studied all family-independent discrete symmetries up to Z 4 . They found that only 
R-parity and baryon-parity are anomaly-free in this sense. Baryon-parity has the 
advantage of forbidding the dangerous dim-5 AB ^ operators (|3|). Both symmetries 
are consistent with all known experiments but lead to drastically different predictions 
for future collider searches of supersymmetry 0, |10| . 

1 We have used the same symbol L for lepton number and for the left-handed lepton doublet. The 
context should always make it clear which is meant. 
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Global symmetries suffer from the same gravitational breaking as discrete sym- 
metries. However, they can not be the remnant of an anomaly-free broken gauge 
symmetry. In addition a global baryon number for example leads to problems with 
baryogenesis and if broken a global symmetry can lead to an unwanted axion. There- 
fore we do not further consider global symmetries. 

Thus it is most likely that the proton is protected by a (broken) gauge symmetry. 
A first solution to the problem of baryon-number violation via a gauge symmetry 
G = U(l)x was considered by S. Weinberg |2j. His approach was to give all the matter 
fields charges Qx of the same sign and the Higgs fields the opposite sign.0 In order to 
guarantee an anomaly-free theory it is then necessary to include fields which transform 
non-trvially under Gsm, m particular also additional (577(3)) coloured fields. Similar 
models have later been constructed by ||12|| . However, the extra coloured field must be 
massive and such a mass term is in general not U(l)x gauge invariant. In principle 
it can however get a heavy mass when U(l)x is broken. 

Besides the problem of giving Gsm non-singlets a mass such fields also drastically 
affect the success of the unification of the gauge couplings in supersymmetry provided 
their masses are below the unification scale. We shall thus only consider G^M-singlets 
as additional fields to cancel the U(l)x anomalies. In this case there is no longer an 
anomaly-free U(l)x where all matter fields have the same charge and the problem of 
AB,AL 7^ interactions must be reconsidered. 

This approach was first studied by Font et al. | JI3|j . The authors searched for a 
U(l)x which prohibits all dangerous dimension four terms (||]). They found that U(l)x 
necessarily acts on the MSSM matter fields as a linear combination of hypercharge 
Y and the third component of right-handed iso-spin 1%. Since (B — L) is a linear 
combination of Y and 1% , U(l)x can also be considered as a linear combination of Y 
and (B — L) 

X = (XyY + a (B _ L) (B — L) + a s S, (5) 

here S is the part of U(l)x acting on the singlets. Later, in || || this work was 
extended to search for a U (l)x which only prohibits a subclass of the terms (0). They 
found both a gauged baryon-number and a gauged lepton-number. However, the 
anomaly-free solutions require additional fields (beyond those of MSSM) transforming 
non-trivially under Gsm, for example the leptoquark of E 6 . As stated above, we do 
not further consider this approach here. 

Since in supersymmetry the proton is most likely protected by a (broken) gauge 
symmetry, we study the possible U(l)x solutions in detail. We extend the work of 
in several points. First, we determine the most general U(l)s of ([5]) which is anomaly- 
free. For three additional Gsm singlets we find an infinite set of rational solutions. 
The singlet charges are generally not identical and thus an interpretation in terms of 
right-handed neutrinos and an embedding in (a family independent) SO (10) is not 
necessary. Second, due to the importence of gravitational effects when considering 
discrete symmetries we embed our models in local supersymmetry. We consider local 
supersymmetry as an effective theory of a more complete unified theory including 

2 This is just the gauged Fayet U(l) which was introduced in order to give the scalar fermions 
large positive mass corrections after supersymmetry breaking. 
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gravity. We can then study the spontaneous breaking of supersymmetry and U(l)x 
explicitely. We determine the scale of U(l)x breaking (Mx ) dynamically in terms 
of the Planck scale Mp\ and the supersymmetry scale. Before U(l)x is broken, but 
below Mpi the superpotential can in principle contain all non-renormalizable terms 
compatible with Gsm <8> U(l)x suppressed by powers of Mp]. These terms could 
be generated by loop effects from the broken unified theory. We consider all such 
terms and their possible effects at low energies, eventhough in any explicit (yet to be 
constructed and therefore unknown) unified model we expect only a subset of these 
terms to be generated via loop effects. After the breaking of U(l)x at M x , the 
renormalizable superpotential will then always contain all terms (0) suppressed by 
powers of (M x /M P i). We thus predict the scale of B and L violating effects including 
neutrino masses in the effective theory and find them compatible with experiment. 
An interesting effect of this scenario is that we always find an unstable but long-lived 
lightest supersymmetric particle (LSP). We find the cosmological constraints on the 
LSP lifetime to be significantly stricter than those due to proton decay and only barely 
compatible with our models. 

The above discussion refers to U(l) gauge symmetries. In a separate publication 
we discuss gauged U(l) R-symmetries [14]]. 

2 Anomaly Cancellation Conditions 

We shall embed our Gsm <2> U(l)x models in N = 1 local supersymmetry. In accor- 
dance with our philosophy stated above, the matter chiral multiplets are those of the 
MSSM with the addition of Gsm singlets, N, and Z m . These multiplets are denoted 
by 

L: (1,2, -i,/), E : (l,l,l,e), Q: (3,2, 

U: (3,1, -§,«), D: (3,l,±,d), H : (l,2,-\,h), (6) 

H: (1,2, 1 h) N : (1,1, 0,n), Z m : (1,1, (W), 

where we have indicated in parentheses the Gsm <8> U(l)x quantum numbers. We 
shall assume that the superpotential in the observable sector has the form 

g (0) = h%UE cj H + h%Q l D cj H + h%C?U cj H + h N NHH (7) 

where He, ho, hu h<N and are Yukawa couplings. Thus at this stage we assume the 
theory conserves R-parity. We have added the term NHH instead of fiHH as in 
the MSSM, in order to incorporate a possible solution to the /z-problem. We shall 
show below that this is not possible. The singlet couplings will be determined by the 
charges z m obtained from the solutions to the anomaly equations below. 

To build a realistic model the superpotential should be gauge-invariant and the 
new gauge symmetry U(l)x should be anomaly-free, i.e. 

l + e + h = 0, (8) 
q + d + h = 0, (9) 
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= o, 


(13) 


3[2l 3 + e 3 + 6q 3 + 3m 3 + 3d 3 ] + 2h 3 + 2h 3 + n 3 + £ z 3 m 


= o, 


(14) 


A A A 


= o, 


(15) 




= o, 


(16) 


3[2Z + e + 6g + 3m + 3d] + 2{h + h) + n + ^ 2 m 


= 0. 


(17) 



The last six equations are the Y 2 X,YX 2 ,X 3 , (SU(2) L ) 2 X, (SU(3) C ) 2 X and gravita- 
tional anomaly |y| equations, respectively. We have assumed that the U(l)x -charges 
are family independent, e.g. h = h = h = I an d there are three generations. In solv- 
ing this set of ten equations we first notice that the Eqs.(|8])- ([T0|) , (|12| , [1~3D , (0), and 
( ITS ) are independent of the singlets N, Z m . The solution to these seven equation can 
be expressed in terms of two variables which we choose to be I and e 

12 4 
h = —I — e, h = I + e, q — — I, u = — I — e, d — -I + e. (18) 

3 3 3 

Inserting the values for h and h into Eq.(|ll|) we obtain 

n = 0. (19) 

Since the U(l)x -number of the singlet N is zero it does not affect the anomalies 
and can be discounted. We thus eliminate it from our further discussion and replace 
NHH in #(°) by 

fiHH. (20) 

We thus offer no new insights on the /i-problem. However, there is also no Peccei- 
Quinn axion in the theory. The remaining equations involving the singlet field charges 
are 

3(2Z + e)+]Tz m = 0, (21) 
3{2l + ef + Y.zl = 0. (22) 

The choice I = — e/2 would give the chiral multiplets in the observable sector U(l)x 
-numbers which are identical to the Y-numbers. Thus, non trivial solutions (U(l)x ^ 
U(1)y) are only possible through the inclusion of the singlets Z m . We also can see 
from the above equations that with the minimal field content (no singlets) charge 
quantization (q d = q e /3 etc.) is given as the unique solution to the anomaly equations. 
This is true in the SM and MSSM and independently of the gravitational anomaly 
equation |[6L |TJ]. 



The cases of I = — e/2 with non-trvial singlet charges are not very interesting and 
we shall require 21 + e 7^ 0. For only one singlet there is no solution to Eqs.(ETip2|). For 
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two singlets there are no real solutions. The only real solution of the cubic equation 
in fl22|) has e = —21. For three singlets an obvious solution is 



zi = z 2 = z 3 = -(2/ + e). (23) 
We shall normalize Z\ — l since the equations contain no constant term. For 

Z2 = ^(x + y), z 3 =~(y-x), (24) 
we then obtain the quadratic equation in x: 

21 + e = ~(y + l), (25) 
27yx 2 + (y-2) 2 (5y + 8) = 0. (26) 

There are two classes of solutions. For y = 2 we must have x = 0; this corresponds 
to the charges 

21 + e = -1, Zl = Z2 = Z3= 1, (27) 

This is the same as (^3|). For I = — 1 this corresponds to the gauged (B — L) symmetry 
of SO (10) and the three singlets are interpreted as the right-handed neutrinos of the 
1650(10)- 

The other class has y ^ 2. The condition for a real solution of Eq. (|26|) is ye [— ~, 0). 
In order to have a rational solution we must have 



5y + 8 2 



3y ° r y = -5W qeQ - (28) 

Then 

x = ±*(— + 2). (29) 

The two signs for x correspond to the interchange z 2 z%. Choosing the plus sign 
and solving for the singlet charges we obtain 

(»-!)(«■ + « + 4) (» + l)( g '-» + 4) 

' 5 + 3g 2 ' 3 5 + 3g 2 V ; 

The charges for the standard model fields are then determined by two free parameters 
/ and q via 

1 -q 2 

e = ^7 - 21. (31) 

5 + 3q 2 K J 

We have thus obtained the complete set of anomaly-free solutions for U(l)x and three 
additional singlets.0 When discussing the details of models we shall mostly focus on 
the specific solutions 

1 4 

q = 0, 2l + e = -, zi = l,z 2 = z 3 = --, (32) 

1 5 

q = 3, 21 + e = --, z x = z 2 = 1, z 3 = --. (33) 



3 The case 21 + e = is included as q — ±1, z\ = 1, z% = 0, Z3 = — 1. 



6 



The other models lead to very similar results as we shall see. We can use these 
solutions to determine the coefficients a» of (0). The additional ?7(1) 's which only act 
on the singlets are given by 

U(l) Sl : (z u z 2 ,z 3 ) = (1,-1,0), (34) 

U(1) S2 : ( Zl ,z 2 ,z 3 ) = (1,0,-1), (35) 

and we have 

U(l) x = 2(! + e)U(l) Y -(2l + e)U(l) {B - L) (36) 

+ ^(l-2z 2 + z 3 )U(l) Sl + ^(z 2 -2z 3 + l)U(l) S2 . (37) 

This contains all solutions, including fl2T|). For the solutions ( |5CD the a, are given in 
terms of /, q 



v 5 + 3g 2 y (D ^> 5 + 3g 2 ' 
g 3 - g 2 + 3g - 3 g 3 + q 2 + 3g + 3 
a * = sT^ 2 " ' aS > = 5 + 3g 2 • (38) 

Thus we have found the following result. Due to the free parameter / above, and taking 
the Z \ = 1 normalization into account, any rational linear combination of U(l)y and 
U(1)b-l can be made anomaly-free acting on the MSSM fields with three additional 
Gsm singlets. For each linear combination of U(1) Y and U(1)b-l there is an infinite 
set of charges zi,z 2 ,zs given by ( p0|) which lead to an anomaly-free Gsm ® U(l)x 
model. These can be expressed as specific rational linear combination of U(l)s lt s 2 - 
Next we dynamically break supersymmetry and U(l)x- 



3 Breaking of Supersymmetry and U(l)x 

To have a realistic model both supersymmetry and U(l)x must be broken at low 
energies. Since we have a locally supersymmetric theory, it is possible to break super- 
symmetry spontaneously. The easiest way is to utilize a hidden sector whose fields are 
singlets with respect to the Standard Model gauge group. Depending on whether the 
U(l)x -symmetry and supersymmetry are to be broken simultaneously or not, these 
singlets would have or not have non-trivial U(l)x -numbers. 

We can break supersymmetry independently of the U(l)x by adding to the system 
one singlet Z with zero U(l)x -number. This clearly has no affect on the anomaly 
equations. Then we can take the superpotential 

g = m\Z + (3) + g\Z 1 , Z 2 , Z 3 ) + g^fa) + g"{Z h Sj), (39) 

where g(°'(Si) is the observable sector superpotential which only depends on the SM 
chiral superfields S^. g"{Z^Sj) is a non-renormalizable part of the superpotential 
involving also interactions between Z^ and Sj. We will discuss it in more detail in the 
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next section. The first term is the Polonyi potential, where j3 is a constant. It will be 
fine-tuned so that the cosmological constant is zero. 

In the following we shall study the two models fl3"2] , |3"B| ) from the previous section. 
The other models will give similar results as long as two singlets have different sign 
charges. This is always the case for (|30"1). It is not the case for the solutions (|23|) 
which correspond to U(l)x = —(2l + e)(B — L). In this case it is not possible to write 
a superpotential of weight zero for the singlet fields. 

For the model ( |33"D it is not possible to write a renormalizable potential as a 
function of Zj, Z 2 , Z 3 . The simplest function we can write is 



g> = k% 4 ( £ a p Zf P Z p 3 , 

\P=0 



(40) 



and we can assume the symmetry Z% <-> Z 2 which implies that a p = a§_ p . We 
have inserted the Planck scale in (|40| ) in order to avoid introducing new scales, k = 
\/8iyGn ~ 1CT 19 GeV~ l . For the model ([32]) the superpotential is given by 



g' = K'Zt a p Zf P Zl 

p=0 



(41) 



The two models lead to very similar results. 

We now analyze the potential V(Z, Z±, Z 2 , Z 3 ) with g' given in Eq.(^). We shall 
take the Kahler function to correspond to minimal kinetic energy for the scalar fields 



/C 



K 



{\Z\ 2 + \Z t \ 2 + \Z 2 \ 2 + \Z 3 \ 2 + \Zi 



(42) 



If we ignore the low-energy fields Si, i.e. consider g", g(°'<^. < Z > 3 , < Zi > 3 we obtain 
for the potential 



V 



+ 



+ 



+ 



rn 



+ —Z\m 2 (Z + (3) + n 2 Zl ]T a v Z\~ v Z, 
^ ' p=0 



K 



K 



■<Zl £(5 - p)a p Zf P Zl + £-Zi \m 2 (Z + f3) + k*Z\ £ a^Zfi 
P =o V P =o J 



<Zt J2 pa p Zt P Z p 2 ~' + y ^ 2 ( m 2 (Z + /3) + k"Z\ £ a^Zf 



P =l 



K 2 - 



Ak 6 ZI a p Zl- p Z p 2 + —Z 3 I m 2 (Z + /3) + K e Z^ ]T a p Z\~ v Z? 



p=0 



p=0 



— —K 



m 



[Z + {3) + n'Zt £ a p Z\- p Z p 2 

p=0 



,9 ■ 



\Zi\ 



(43) 
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g' is the U(l)x gauge coupling, Z is the complex conjugate field of Z. We have 
assumed here that the gauge field kinetic energy term is minimal: f a p = 8 a p. Since the 
fields Zi transform under the extra U{l)x group the corresponding D-term appears 
as the last term in the above potential. There are many possible minima for this 
potential, mainly with < Z\ >~< Z 2 >~< Z 3 >m< Z >m A. However, we are mainly 
interested in the situation where <Z>^ - ~^><Zi>, <Z 2 >, <Z 3 > and we shall tune 
(3 and the constants a v so that V is zero at such a minimum and positive definite, and 
where the condition 



|Z 1 | 2 +|Z 2 | 2 ~-|Z 3 | 2 (44) 

is satisfied. In this case we can expand around the vev Z m - and write down the 
effective potantial as a function of Z±, Z 2 , and Z 3 



V 



+ ml /2 \Zi\ 2 + m 3/2 [Z l g,i + (A- 3)g + h.c. 
2 5 2 



(45) 



where g is the same as g' up to a multiplication factor which can be absorbed in the 
a p . The parameters A and 777.3/2 are related to the potential g: m 3 / 2 ~ k 2 <g>. For 
the Polonyi type potential 7713/2 ~ Km 2 . It is now easy to minimize the potential (f45|). 
For simplicity we shall assume that <Z 2 >= 0. Then 



V 



K 



12 



lQlZ^lZA 10 + 25\Z 3 \ ti \Z 1 \ 



+ 777 



+ K 6 A'm 3/2 \Z°Zl + h.c 
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Zi\ 



3/2 



4' 



\Zi\ 2 + \Z 3 \ 

■ |2 



minimizing with respect to Z\ and Z 3 one finds that the equations imply 



(46) 



\Zi\ 

and < Z\ > satisfies the equation 

576 J 



0(m 2 3/2 ) 



x 2 + 



j m 3/2 x + - ^ 



32 



m 



3/2 







where x = k 6 < Z\ > 7 . The solution of this quadratic equation gives 

' "^3/2 



o 



100 K< 

10 16 GeV 



M 



x 



(47) 



(48) 



(49) 
(50) 



From the above analysis we deduce that for the solution with Z\, Z 2 symmetric this 
gives 

5, 



<Z X > 2 =<Z 2 > 2 



\Z 3 \ 2 + 0(m£ /2 ) 
0(M 2 X ) 



(51) 
(52) 
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so the mass for the U(l) gauge bosons is O(Mx)- The approximation we used to obtain 
this solution is valid because shifts in the low-energy sector due to the breaking of the 
U(l)x symmetry are of order 

n 2 g' « 0{k 8 M 9 x ) = 0(1 GeV) (53) 

For the solution ( [23] ) there is no superpotential and the above mechanism will not 
apply. It is still possible to break supersymmetry and U(l)x by taking a non-minimal 
kinetic energy of the form 

3 3 

EW + al^lVhEW (54) 

i=l i=l 

However, this would only lead to solutions < Z\ >=< Z 2 >=< Z 3 >= O(^). As will 
be clear from the discussion in the following section these solutions will lead to terms 
in the effective potential which give rise to too fast proton decay since Mxk = 0(1). 
It does not seem that one can construct a realistic model in this case. 

Since the U(l)x -symmetry is broken, there will be U(l)x -breaking terms in the 
effective action, and these will be induced through the U(l)x -gauge and gaugino 
interactions. Such terms can lead to proton decay, or to baryon and lepton number 
non- conserving processes. We discuss them in the next section. 



4 Baryon- and Lepton-Number Violation 

The baryon- and lepton-number violating terms of Eqs. (fJH) and can lead to proton 
decay. We first discuss the dimension four terms (fj). Given our solutions from Section 
^| their charges Qx can be expressed in terms of I, e, 

Qx(LLE, LQD,UDD, LH) = 21 + e. (55) 

For a non-trivial U(l)x we had required 21 + e ^ and thus at tree-level all dimen- 
sion four terms are excluded. However, local supersymmetry is a non-renormalizable 
theory. We thus expect our model to be the effective action of a more complete uni- 
fied model at a higher scale. We expect this unified theory to include AB, AL ^ 
effects. The non-renormalizable terms in our model below the Planck scale will then 
be obtained as the 1-loop effective action of this unified theory. However, for lack of 
knowledge of the unified model we shall consider the most general non-renormalizable 
interactions between the singlets Zi and the MSSM fields Si which are Gsm <£> U(l)x 
gauge-invariant 

3 

K N Z n 1Z n 2Z nsJJ S , (gg) 
i=l 

where N = ^1+^2+^3 and the Si product is a trilinear term of Eq.(^|).0 After breaking 
U(l)x in principle one has to find the one-loop effective action below M, which will 

4 Since the low-energy effective Lagrangian always contains the term HH, and the terms 
LH, LLE, LQD are either present or absent simultaneously, it is possible in our models to rotate 
away the LH terms and we do not further consider them. However, we emphasize that for other 
models this must always be checked in each case. 



10 



include AB, AL ^ terms. As our singlet interactions are non-renormalizable from 
the start, a calculation of the effective interaction is not very reliable. But all terms in 
the effective action are necessarily derivable from non-renormalizable terms of the form 
(|56|). We shall then consider the terms obtained from (^) by replacing Z^ — >< Zi>. 
If all the AB, AL ^ effects generated in this manner are suppressed, then we can 
conclude that AB, AL ^ effects are not observable in the laboratory. 

When replacing Zi —>■< Zi> we obtain an effective superpotential containing the 
terms 

(KM x ) N Y[Si (57) 

%=\ 

with effective coupling constants = (nM x ) N < 1. Here M x ~< Z { >^ 10 16 GeV 
is the scale of U(l)x breaking. Since the terms LQD and UDD have the same charge 
Qx they will be suppressed by the same power of A ~ e . If we have an effective 
superpotential 

g e ff(Si) = XLLE + X'LQD + \"UDD (58) 

then we can use previous proton decay rate calculations and the experimental 
lower bound on the proton lifetime jl9) to estimate the extremely strict bound 

A . A , <ur ,Q^. (59) 

Since e ~ we must have iV > 5. In our first model (21,22,23) = (1, — I, - 5), 
21 + e = I and the lowest dimensional term is 

K 7 Zf(Z 2 Z 3 ) 2 g e "(S t ), (60) 



and indeed we get the appropriate suppression with N = 7. For our second model, 
(21, 2 2 , 23) = (1, 1, —4), 21 + e = — \ and the lowest dimensional term is 

K 7 ZlZt g^iSi), (61) 

and again N = 7. The general case for rational q is not soluable. However, we can 
determine an upper bound on the suppression. For terms symmetric in Z2, Z% we 
obtain 

(Z,Z 2 rzr g eff (S t ) (62) 



K 



2ni+n 3 



Q = ^r 2 +^ + ^- 2 ^o (63) 



which has the charge 

5 + 3g 2 3 5 + 3g 2 
and must vanish. For q 7^ 1 the lowest dimensional solutions are n\ — 3, 713 = 1 and 
rii = 2, n 3 = 3, and thus N = 7 once again. However, it is possible that for certain 
values of q the suppression might be weaker. 

We now turn to the dimension five terms. Their charges are given by 

Qx(LHHH, HHE*) = ^Qx(LLHH) = 21 + e (64) 

QxiQQQH, HQtJE, H*, QtJL*, UD*E) = -(2/ + e) (65) 
Qx(QQQL,UUDE) = (66) 
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For the terms in Eq(Q,^) the discussion is identical to that of the dimension four 
terms above. They are suppressed by e 7 in any effective theory. The last two terms 
are B — L invariant and they can not be excluded by any gauged U(l)x, given our 
assumptions. It is however not so clear whether these terms actually pose a problem. 
In Ref. f2j the partial proton decay rate via the operator (QQQL) F was estimated to 
be 

r pW Q Qi) «p) 2 g^. (67) 

Here M is the mass suppression of the non-renormalizable dimension-5 term in the 
effective Lagrangian. For us M = e is the electric coupling which enters when 
the final state s-fermions are converted via electroweak gaugino exchange to their R- 
parity even partners. This involves a 1-loop diagram, hence the 87r 2 . fh is an effective 
supersymmetric mass, e.g. rh = ™ g augino/ m2 squa rk- ^ i s the coupling constant of the 
effective operator; it is squared since when generated from renormalizable interactions 
it involves at least two coupling constants. The dominant decay mode corresponding 
to (^) is p — > K + u ||; the experimental bound is []TS]] r proton > 10 32 a. This then 
corresponds to 

m n " 2 



* <10 IToog^J ' (68) 

This is the same order as the muon Yukawa coupling in the Standard Model, which 
is not explained but is generally also not considered to be unnatural. None the less, 
we proceed to consider the conditions where these terms are suppressed. In order for 
the terms (|66|) to be induced at the loop level when U(l)x is broken, they must occur 
through the exchange of Z$ superfields. This implies these terms are of the form 

f(Z 1 ,Z 2 ,Z 3 )(aQQQL + (3UUDE) (69) 

a,P are constants and f(Z 1 ,Z 2 ,Z 3 ) must have U{l)x charge zero. The lowest order 
contribution is for f(Zi, Z 2 , Z 3 ) to be of the same form as g' in ( flCf ) 

f(Z u Z 2 , Z 3 ) = k 10 ZI £ Z\- p Zl (70) 

p=0 

and this gives rise to the effective coupling 

G 5 =< / >~ K 10 < Zi > 9 ~ I0^ 36 GeV' 1 (71) 

which is highly suppressed. A constant term corresponds to a tree-level term and 
is not allowed. Thus the dimension five terms are safe as well. However, when the 
unknown unified theory is broken, and before U(l)x is broken the QQQL term could 
be generated directly in the effective theory. In this case the coupling must satisfy 
the above bound (1681). 
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5 LSP Lifetime and Neutrino Masses 



We now turn to two applications of our analysis. We have found that in general we 
expect an effective i?-parity breaking superpotential fl5§D with highly suppressed but 
non-zero couplings 

A « A' « A" « 10~ 21 . (72) 

The terms of g e ^(Si) induce the decay of the lightest supersymmtric particle (LSP).0 
For a photino LSP the decay rate has been calculated by Dawson PU] 

r- - ^LlA (73) 

a is the fine structure constant, A is the coupling in g ef f and Mj is for example the 
scalar electron mass. The more general decay of a neutralino LSP is given in flQ |. 
The lifetime of the LSP is then given in natural units as 



23 / 10- 21 \ (50GeV\ ( Mf 



T -- = w a [^r) {—) [wg#) • < 74 > 

The large difference in lifetime compared to the proton is due to the different masses 
and the two powers less dependence on A 

/ \ 5 2 



Tproton \ ^'proton J V A 



10 49 (75) 



for = 50 GeV; e is the electric charge quantum. 

If we allow a supersymmetric mass range up to 1 TeV for the scalar fermion mass 
and LSP masses up to Mj we obtain in terms of the lifetime of the universe t u = 
2 • 10 17 s 

10 3 < H < 10 7 . (76) 

There are several large uncertainties in this result. First, we have taken the singlet 
vev < Z >= O(10 16 GeV). The insertion of a simple factor of three either way changes 
A = e 7 by a factor 10 ±3 and by a factor 10 ±6 ! Second, we have only considered 



a photino LSP. As discussed in Ref.|T0| in detail the LSP lifetime can vary by many 
orders of magnitude for a general neutralino LSP, depending on the MSSM parameters. 

Surprisingly enough there are severe constraints on a long-lived but unstable LSP 
pl| , |2~2|| , even for > t u . The relic density Q^h 2 of a stable LSP is typically in 
the range — 10 [23|; if constrained to the minimal N = 1 local supersymmetric 



model it is 10~ 2 — 10 |[24|| . For as large as 10 7 r u only one in 10 7 LSP will have 
decayed today. However, for such large relic densities this is still a large number 
and the decay products can lead to observable effects. In particular final-state decay 
neutrinos can be observed in the large underground detectors. For a LSP mass below 



5 The allowed dimension 5 terms of (JfS^) conserve R-parity and are thus irrelevant for the decay 
of the LSP. 
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1 TeV the strictest bound [211] comes from the experimental upper limit on upward- 



going muons from the 1MB detector |25j . In terms of the LSP branching ratio B v to 



muon neutrinos this bound can be expressed as [2~I| 

^ < 4 ■ lO- 9 ^- 1 ! ln-^l + I ^ 7 ) (77) 

Assuming all operators in fl5"B|) have equal strength the muon-neutrino branching 
fraction is of order 1/5. The lifetime is reduced by about a factor of 35, assuming 
< m top .[] For an LSP mass of 150 GeV we thus obtain 

n,h 2 < 6 • 10~ 9 ^, (78) 



and we see that only the upper range in Eq.flTq) is still allowed. This analysis is a 
bit crude in many respects. For example, both the relic denisty and the LSP lifetime 
depend strongly on the MSSM parameters and it would be appropriate to perform a 
correlated analysis. This is beyond the scope of this paper. 

Thus although the LSP lifetime is many orders of magnitude shorter than the 
experimental bound on the proton lifetime ij proton > 10 22 r n ) the resulting bound on 
the couplings A, A', A" of ( |58D is actually 2 — 3 orders of magnitude stricter, if all 
operators are of equal strength. 



The non-renormalizable interactions ( pB[ ) can also give rise to neutrino masses. 
The singlets Zi can couple to the left-handed neutrinos via 

k 6 LH Y[ Zi (79) 

i 

leading to Dirac neutrino masses of order e 6 <H>. The leading Major ana neutrino 
mass for the singlets is given in (fermionic) analogy to Eq.(^); it is of order 7723/2 = 
(9(100 GeV). We then obtain the see-saw mass matrix 



e 6 <H> 
e 6 <H> m 3 / 2 



$0) 



The light-neutrino mass is given by 

(V 6 <H>) 2 

m v w 1 J w 3 • 10~ 18 eV, (81) 

fn 3/2 

which is extremely small and unobservable. This is consistent with all experiments 
but can not explain the hints for neutrino masses such as the solar neutrino puzzle. 



6 Conclusion 

We have found the most general anomaly-free U(l)x with three Gsm singlets. In 
local supersymmetry the breaking of U(l)x is predicted to be at the scale 10 16 GeV. 

6 For m-y > mtop these numbers change to 1/6 and 45; the product only changes from 7 to 7.5. 
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The most general gauge-invariant non-renormalizable superpotential is consistent with 
proton decay experiments. It is only barely consistent with the cosmological bounds 
on the decay of the lightest supersymmetric particle. This poses a severe constraint 
for all supersymmetric models. 
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